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A TURBULENCE INDICATOR UTILIZING THE DIFFUSION OF HEAT 



By G. B. SCHUBAUBB 



SUMMARY 



In cooperation with the National Advisory Committee 
for Aeronautics, the National Bureau of Standards for 
several years has been studying methods of determining 
the turbulence in wind tunnels, especially by the "hot-wire " 
method and the sphere method. This paper describes a 
third method. 

The effect of turbulence upon the diffusion of heat from 
a small electrically heated wire in an air stream was 
investigated. The turbulence of the stream was introduced 
by a series of geometrically similar screens placed one at 
a time across the upstream section of the tunnel. With the 
wire set at various distances from the screens, curves of 
temperature distribution were obtained by traversing the 
heated wake at a distance of 2 inches behind the wire with 
a small thermocouple. A single relation was found to 
exist between the width of the wake at half maximum and 
distance in screen wire diameters from the several screens. 
The correlation of width at half maximum with percentage 
turbulence, as measured by the "hot-wire" method, could 
be represented approximately by a single curve. 

INTRODUCTION 

The evaluation of wind-tunnel turbulence is by no 
means common practice in aerodynamic laboratories, 
even though turbulence is recognized as an important 
aerodynamic factor in wind-tunnel tests. This situa- 
tion is due largely to the fact that a turbulence deter- 
mination is attended with more difficulty and uncer- 
tainty than is the average aerodynamic measurement. 

It is possible to determine quantitatively the ratio of 
the root-mean-square of the speed fluctuation at a 
point to the average speed. This quotient times 
100 is usually called the percentage turbulence. 
The percentage turbulence may be measured by a mod- 
ified hot-wire anemometer and accessory equipment, 
including an amplifier and an electrical network to 
compensate for the lag of the wire (references 1 and 2). 
The apparatus is expensive and cumbersome and its 
operation requires considerable skill and care on the 
part of the operator. 

Simpler turbulence indicators, which, if desired, might 
be calibrated in terms of percentage turbulence by 
comparison with a "hot-wire" instrument, have been 
sought. A device of this sort, which has found wide 
use in comparing the turbulence of different wind tun- 



nels, is the sphere (references 3, 4, and 5). The corre- 
lation of its critical Reynolds Number (defined as the 
Reynolds Number for which the drag coefficient is 0.3) 
with percentage turbulence has been made in several 
laboratories, but the universality of this correlation is 
still. under investigation. 

The present article describes a scheme of turbulence 
measurement based on the diffusion of heat from a 




Fiarmi 1.— Diagram or apparatus used to traverse the heated wake. 

heated wire. It has been found by experiment that 
the width of the heated wake increases as the turbu- 
lence increases, indicating a relation between the rate 
of diffusion and the degree of turbulence present in the 
air stream. 

APPARATUS 

The apparatus for producing the heated wake and 
traversing it to determine its width is shown in figure 
1. A platinum-iridium wire 2.8 inches long and 0.002 
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inch in diameter is located at A. At B two no. 36 
wires, one of copper and the other of constantan, are 
soldered together forming one junction of a thermo- 
couple. The other junction of the same size at C is 
placed just far enough from B to be completely outside 
the heated wake. The prongs carrying the heating 
wire and the thermocouple are attached rigidly to the 
arm D, which can be rotated about a vertical axis 
passing through the heating wire to carry the thermo- 
couple B across the wake. The platinum-iridium 
wire was heated by an electric current ranging from 
0.2 to 0.8 ampere, depending on the heating desired. 
The current could be maintained sufficiently constant 
by means of a rheostat and a 0 to 1 ampere range 
ammeter. The thermocouple was connected to a 
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Fioube 2.— Distribution carves obtained by traversing the heated wake at a dis- 
tance of 2 Inches from heated wire. Cnrve (a): LI percent turbulence. Curve 
(b): 3.5 percent turbulence. 

sensitive wall galvanometer (sensitivity 0.06 micro- 
volt per mm per m scale distance). 

Figure 2 shows two typical distribution curves 
obtained with the apparatus for a spacing of 2 inches 
between the heating wire and the thermocouple. 
Curve "a" corresponds to a turbulence of 1.1 percent 
and "b", to a turbulence of 3.5 percent. For both 
curves the heat input was 5.7 watts and the air speed 
36 feet per second. Since the maximum ordinate of 
curve "a" corresponds to a temperature rise of only 
3° C, and the copper-constantan thermocouple has a 
thermoelectric power of 41 microvolts per degree 
Centigrade, the necessity for a rather sensitive gal- 
vanometer is obvious. 

THEORY 

It is desirable that the heating wire have a small 
diameter so as to create little disturbance of the flow. 
The transverse flow of heat will be due to ordinary 



molecular conduction plus the eddy conduction arising 
from the turbulence initially present in the air stream. 
Let us consider first the case where the stream is com- 
pletely free from turbulence and the wire produces a 
negligible disturbance. Regarding the wire as a line 
source of heat and the air speed as uniform, the tem- 
perature distribution is given at points not too close to 
the wire by the equation: 1 

AT=AT mar e~^ (1) 

where y is distance across the wake measured from the 
center; x is distance downstream from the wire; AT 
is the temperature rise at any point (x, y) ; A T max is 
the maximum temperature rise at the point (a;, o); v is 
air speed and p, c, and k are density, specific heat at 
constant pressure, and thermal conductivity of the air, 
respectively. 

It will be observed that equation (1) is anal- 
ogous to the expression for the linear flow of heat 
from an instantaneous plane source. The approxi- 
mations involved when applied to the present prob- 
lem may be illustrated by visualizing the process 
described by equation (1). A line of air particles, after 
being heated by contact with the wire, is carried on 
in the x direction with the speed v. While the line 
moves through the distance x, during the time interval 

— > the heat flows outward from the line of particles in 

the ±y direction, and at x has attained the distribution 
given by equation (1). Clearly, the finite size of the 
wire, the conduction of heat in the x direction, the 
retardation of the flow by the wire and the stirring 
action due to whatever eddies are formed in the wake 
of the wire have been neglected. 

If the air incident upon the wire is turbulent, there 
is an apparent conductivity expressed by (k+P), 
where /S is the eddy conductivity due to the turbulence. 
Equation (1) then becomes 

AT=AT max ^ k +^ 1 (2) 

The width of the wake at half maximum tempera- 
ture is selected to represent the width characteristic 
of the wake. From curves such as shown in figure 2 
the angle subtended at the wire by the wake width at 
half maximum may be obtained. For convenience, 
the results will be expressed in terms of this angle 
rather than in actual width, the term angular width 
and the symbol a being used to designate the angle. 

' The exact form of the temperature distribution behind a line source In lamlnnr 
air flow as given by Drew (Trans. Am. Inst, Ohem. Engineers, vol. 28, 1931, p. 30) la 

whjre q Is the heat output P3r unit length of the wire, i stands for the combination 
and -Koto is a Bessel function of xero order. At large values of t, 

Ko (i) may be set equal to -^/^«-*. After expanding -/r»+F" for >«* and neglect- 
ing small terms, equation (1) follows. 
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When ^ — =?r in equation (1), y is the half width 

OEqS 

57.3 



at half maximum. Making the substitution 2y=jr^ 



and solving for oo, we obtain 
ao=190.8-^ 



k_ 

p cvx 



(3) 



where a 0 is the angular width in degrees at half maxi- 
mum for laminar flow. Making a similar substitution 
in equation (2), we obtain 

I Vc+0) 
p cvx 



a= 190.8-^ 



(4) 



RESULTS 

Measurements were made in the 4^-foot wind tunnel 
of the National Bureau of Standards, in an air stream 
made turbulent by a series of geometrically similar 
square-mesh screens placed one at a time across the 
upstream section of the tunnel. The amount of tur- 
bulence obtained from any one screen was varied by 
working at several distances from the screen. The 
first, second, and third columns of table I give, respec- 
tively, the mesh and wire diameter of the screens and 
the distance in wire diameters from the screens for the 
several working positions. The distance was measured 
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Fioube 3.— Variation of a,„v with distance from screens expressed in screen wire diameters. 
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where a[is the angular width for turbulent flow, 
squaring equation (4) 
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p cvx 



p cvx 



= o 0 2 -(-a 2 (uri> 



(5) 



where a, UT t, is the angular width due to the turbulence 
alone. 

From the foregoing expressions, it is clear that the 
heat output does not have to be known.* It is nec- 
essary only that the heating be constant, during a 
traverse of the wake, and of sufficient magnitude to 
allow the temperature distribution to be measured. 

» Tho temperature rise Is so small that tte variation of k with wake temperature 
may te neglected. 



from the screen to a point midway between tho heating 
wire and thermocouple. Since the screen pattern in 
the speed distribution persisted to about 65 wire diam- 
eters downstream from the screen, the measurements 
were made at distances greater than 65 diameters. 

The fourth column of table I gives the observed 
angular widths for a 2-inch spacing between the heat- 
ing wire and thermocouple. The values of a tUTi cal- 
culated by equation (5) are given in the fifth column 
of table I. A comparison of these two columns shows 
that the effect of a<, becomes less as a increases. From 
rims made at wind speeds ranging from 8 to 55 feet 
per second, it was found that a varied with the speed, 
but that a tU rb was constant to within the limits of 



experimental uncertainties. Since a turl> -- 
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the eddy conductivity /S must therefore be proportional 
to the air-stream velocity v. 



cant that the points fall about a single curve in 
figure 3, whereas they follow a family of curves in 



or 
u 

I 

i 





















































( 














I 






















































— H 






\ 


V- 






























t 


k 


5" 






























































X 


\* 


-=/ 


























) 

































































































200 400 600 800 IOOO 1200 1400 
D/d - Distance from screen in wire diameters 



1600 



Fiquee 4.— Variation of percentage turbulence with distance from screens expressed in screen wire diameters. 



Figure 3 shows the variation of a turt with distance, 
expressed in screen diameters, from the several 
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Fioubb 6.— Correlation of o,„ k with percentage turbulence, AT denotes the screen 
mesh. (See Illustration In figs. 3 and 4.) 

screens. Figure 4 is the corresponding diagram of 
percentage turbulence, measured by the "hot-wire" 
method, for the same series of screens. 3 It is signifi- 

' These data were obtained by W. O. Slock, Jr., andtheauthorln connection with 
another research project of the aerodynamio section. 



figure 4. A possible interpretation of this result is 
that the "hot-wire" method differentiates between the 
scale of the turbulence produced by screens of different 
size, but that the diffusion method does not. It is 
reasonable to assume that the eddy size may bear a 
direct relation to the screen size. If so, it seems pos- 
sible that the wire of about one centimeter length 
employed in the "hot-wire" method may average out 
the speed fluctuations to a progressively greater degree, 
and give a smaller response, as the eddies get progres- 
sively smaller. This subject is now under investiga- 
tion at the National Bureau of Standards. 

Figure 5 shows the correlation of a, uri with percent- 
age turbulence. Obviously, from figures 3 and 4 the 
correlation cannot be represented exactly by a single 
curve. However, the majority of points fall near a 
straight fine passing through the origin, the marked 
deviations entering only for the }£-inch and #-inch 
mesh screens. 

When the effect of spacing between the heating wire 
and thermocouple was investigated, it was found that a 
was essentially the same regardless of the spacing in the 
range from % inch to 6 inches. This means that 
the wake is wedge-shaped, the width at half maximum 
varying linearly with x. From these results 0 was 
calculated and found to be proportional to x. 

It is possible to summarize all the results to date in 
an empirical equation for /3 of the following form: 
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where K is a constant of proportionality and is 
the root-mean-square of the speed fluctuation. The 
constant K is equal to 0.63 p c. 

A study of diffusion such as outlined here offers a 
promising field for future research both in theory and 
experiment. 



National Bureau of Standards, 

Washington, D. C, March 8, 19S5. 
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TABLE I 



Mesh of screen (Inches) 


Screen wire 
diameter 
(inches) 


Distance 
rom screen 

in wire 
diameters 


Angle subtended 
at wire by width 
at half maximum 
for 2-inch spac- 
ing between 
thermocouple 
and heating wire 








a degrees 


atmrh 

degrees 






141 


8 99 


3.83 


8 by 5. — 


1.0 


114 
97 


4ies 

6.03 


4.64 
4.90 






06 


6.92 


6.82 






225 


3.10 


2.89 






183 


3.63 


3.36 


3Hby3X 


.625 


155 


3.81 


3.64 




106 


4.85 


4.72 






64 


7.00 


6.92 






739 


1.70 


L28 


1 by 1 


.192 


610 
349 


L85 
2.15 


1.48 
1.84 






216 


2.92 


2.70 






1,540 


1.75 


L35 






1,640 


1.49 


.988 


HbyM 


.093 


1,063 
727 


1.76 
L98 


LS6 

1.61 






460 


2.06 


1.73 






141 


4.46 


4.32 






2,082 


LB 


1.05 


KbyK 


.047 


1,423 

882 


1.63 
1.80 


1.05 
L41 






276 


2.72 


2.48 



Average wind speed=S6.2 ft. (a.)-'-l,104 cm (s)-> 
Average air temperature=18.3° C. 
Average barometer reading"753 mm Hg 
P=1.22l)X10-» grams (cm)"* 
c=0.240 cal. (gram)-i (° C.)-f 
i=6.60X10-' cal. (s)-> (cm)-» [° O. (cm)-»]-i 
I— 6.08 cm 
<w-L114° 



